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' We discuss effects of magnetic fields on proto-neutron star winds by performing numer- 

I ' ical simulation. We assume that the atmosphere of proto-neutron star has a homogenous 

, magnetic field (ranging from 10^'^ G to 5 x 10^^ G) perpendicular to the radial direction and 

QQ ■ examine the dependence of the three key quantities (dynamical time scale Tdyn, electron 

' fraction Y^, and entropy per baryon s) for the successful r-process on the magnetic field 

\ strength. Our results show that even with a magneter-class field strength, ~ lO^^G, the 

, feature of the wind dynamics varies only little from that of non-magnetic winds, and that 

' the condition for successful r-process is not realized. 

^ ■ 

Q . §1. Introduction 

^ ■ Because of their nuclear charges and energetical disadvantages, elements heav- 

• • , ier than iron cannot be produced by thermonuclear reactions. However, their solar 

\ abundances are much greater than those can be obtained in the nuclear statistical 

^ ' equilibrium. One of the processes that contribute to synthesis of this heavy element 

is the r-process (rapid neutron-capture process). With no Coulomb barrier to over- 
come, the seed nuclei of the r-process capture neutrons easily even at low energies, 
forming characteristic abundance peaks at A 80, 130, and 195. The heaviest el- 
ements, in particular, such as ^^^U, ^^^U, and originated only in the 
r-process. 

In the r-process, the typical time scale of neutron capture must be much shorter 
than that of /?-decay. Thus, the neutron-to-seed ratio is a critical parameter for 
the r-process and must be significantly large (> 100) to pile up enough neutrons to 
produce nuclei up to and beyond the third peak of the abundance. Identifying the 
sites that provide with this large amount of neutrons has been an attractive issue 
but is still debated. 

Recent observations of r-process elements in ultra-metal-poor halo stars show an 
excellent coincidence with the solar r-process abundances except possibility for the 
lighter elements. This suggests that the r-process has occurred from early phase of 
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the galactic evolution. Thus, identifying the r-process site should give us information 
on the galactic evolution, the star formation rate, and the initial mass function. 

The studies of the r-proccss elements are also important because some of them 
can be used as cosmological chronometers. For example, uranium and thorium are 
considered to be useful, since they have long radiative half lives (^^^U;4.468Gyr 
^^^Th;24.08Gyr). The age of the metal-poor stars, the lower limit of the age of 
universe, can be estimated by comparing their observed abundance ratio with their 
theoretical initial production ratio. 

Among a number of the r-process sites proposed so far, the neutrino-heated wind 
from the surface of the proto-neutron star (PNS) seems to be the most promising. 
The neutrino-heated wind is an outflow which is considered to take place at the late 
phase of the delayed explosion scenario. 

Since this wind scenario was proposed, the dynamics and the nucleosynthetic 
outcome have been studied. Woosley et al.^) performed numerical simulations by 
using a spherically symmetric hydrodynamics code. They followed up to '^18s the 
wind evolution after the trigger of the core-collapse of a 20Mq star. They found an 
outflow with significantly high entropy per baryon >400A;f,, where is the Boltz- 
mann constant, and succeeded in producing the r-process solar abundance. However 
an outflow with such a high entropy to bring a successful r-process has not been 
reproduced in other numerical simulations.^-*'^) Qian Sz Woosley^^ has done both 
analytical and numerical examination systematically to investigate the properties of 
the wind to seek for a site of r-process. Their analytic models are based on the ap- 
proximate solutions of spherically symmetric steady outflows. They compared them 
with the numerical calculations. They showed clearly that the wind dynamics is 
characterized by three quantities, that is, the dynamical time scale Tdyn, entropy per 
baryon s, and electron fraction Ye- Although they obtained a good agreement be- 
tween analytical and numerical results, the obtained entropy fell short of the require 
value by a factor of 2-3. 

Cardall &; Fuller^) extended the Newtonian analysis of Ref.^^ to the general 
relativistic one. They showed that the increase of mass-to-radius ratio of the PNS 
due to the stronger general relativistic gravity leads to an increase of s and decrease of 
Tdyn, which is favorable for the r-process. Further studies of the wind dynamics in the 
framework of general relativity were done analytically by Otsuki et al.^^ and Wanajo 
et al.,^-* as well as numerically by Sumiyoshi et al.^^ In their works, it was found there 
is a possibility of successful r-process nucleosynthesis to occur with an entropy not so 
high as ^ 400fcft but with a very short dynamical time scale Tdyn ~ 6ms. However, to 
realize this condition, a very small mass-to-radius ratio {GM/(?R ~ 0.3) is required, 
which can not be met by a model with a typical neutron star (NS) mass (~ I.4M0) 
and radius (~ 10km) (see, however, Terasawa et al.^^^). 

Since there had been eventually no successful model with a typical neutron star 
(I.4M0, R ~ 10km), effects of magnetic field on wind dynamics were investigated 
by Nagataki et al.^^ and Thompson.^*') Assuming a monopole-like configuration, 
Nagataki et al. explored effects of the magnetic field and rotation analytically, but 
could not find a suitable condition for the successful r-process. Thompson, on the 
other hand, investigated effects of dipole magnetic field by assuming that the closed 
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magnetic field works to trap the matter, he insisted that the magnetized wind is 
more efficiently heated than in non-magnetic models and claimed that a model with a 
surface magnetic field ^ 6 x 10^^ G would bring a successful r-process nucleosynthesis. 
Though Thompson claimed that he found a successful r-process wind model, he did 
not solve the dynamics, and it is obvious that a more detailed treatment of the 
dynamics of magnetized wind is needed to obtain a firm conclusion. Thus, in the 
present study, we have done it by detailed numerical simulations. As a first step of 
studying magnetic effect on PNS wind, however, we have considered only a radial 
motion of a thin surface layer with a homogeneous magnetic field perpendicular to 
the radial direction. 

The paper is organized as follows. In §2, we explain the numerical method used 
in the present study. In §3, we show the results of our simulations. The results are 
analyzed in §4. Summary is given in §5. 

§2. Numerical method 

We consider the subsonic outflow of magnetized matter from the PNS surface. 
Although this is a genuinely multi-dimensional problem, we treat only the radial mo- 
tion in this paper and solve one dimensional magnetohydro dynamics (MHD) equa- 
tions as explained below. This will approximate the flow near the equator with 
magnetic fields perpendicular to the radial direction. We have employed the code 
employed by Ref.,^) which is an implicit lagrangian code for general relativistic and 
spherically symmetric hydrodynamics (Yamada^^^). In the following, the space-time 
metric is approximated by 

^^2 ^ e2<A(i,m)^2^^2 _ ^2X(t,m) [^^\m^ - r\t,m){de^ + sin" 9 dcf^') , (2-1) 

and we have solved the conservation equations of mass, momentum and energy to- 
gether wit Einstein equations. See Ref.^^^ for detail. Since our aim is to explore 
the propagation of magnetic wind, we have modified the code to include the time 
evolution of magnetic field and its dynamical feedback. 

The temporal evolution of the magnetic field is derived from the equation 

£^F = 0. (2-2) 

Here F = Ff^i^dx^ ® dx'^ is the electromagnetic field tensor, and £u denotes the 
Lie derivative along the four velocity u". The Greek letter stands for 0—3. The 
above equation is valid whenever the electric fields vanish in the comoving frame 
(Achterberg^^)), since the Lie derivative can be written as 

X^F = d(F • u) + (dF) • u. (2-3) 

The first term in equation H2-3() vanishes from ideal MHD approximation, that is 
F ■ u = F^yu'^dx^ = 0. On the other hand, the second term in equation H2-3() 
vanishes independently of the ideal MHD approximation by virtue of the Maxwell 
equation, dF = 0, which is equivalent to F = dA, where A is the vector potential. 
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Since we are using the Lagrangian coordinates, equation (|2-2|) implies that the 
components of the electronmagnetic tensor remains constant for each fluid element, 
which can be written as 

F^^{ti,m)=F^^{tQ,m). (2-4) 

Here, to and ti are arbitrary times. The components of the electromagnetic fields in 
the local inertial frame is obtained by performing the following coordinate transfor- 
mation at each coordinate point 



{t,m) 



(2-5) 



Here, F^is represents the electromagnetic tensor in the local inertial frame. The 
transformation matrix is denoted as ||^, and its components are 



/ g-</'(i,m)g-l 
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r^^{t, m) 











(2-6) 



Thus, from equations ()2-4() - (|2-6|) . the magnetic fields in the local inertial frame at an 
arbitrary time is determined as 



Briti,m) 



r'^{ti,m) 



Br{to,m) 



r(ti,m) 
r[ti,m) 



(2-7) 

(2-8) 
(2-9) 



where Br, Bg, B^ are radial, polar, azimuthal components of magnetic field, re- 
spectively. By taking to as the initial time, the time evolution of the magnetic field 
is determined by equations ()2-7() - (|2-9() . In the present study, we only consider the 
magnetic fields perpendicular to the radial direction and set B^ = B^ = 0. This 
corresponds, for example, to the uniform, dipole, or toroidal field in equator. 

The general relativistic equation of motion and conservation of energy are mod- 
ified by including dynamical feedback of magnetic field. Adding the contribution of 
the magnetic pressure to the original equation, the equation of motion is given by 



'dt 



h 



dm 



dm 



+ 



m 



(2-10) 



where U is the radial fiuid velocity defined as 



(2-11) 
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and r is the general relativistic Lorentz factor defined as 

r^ = l + U^ (2-12 

r 

and r, /i, (7, m, and Pj/ are the inverse of the baryon mass density, the specific 
enthalpy, the momentum transfer from neutrinos to matter, the gravitational mass 
inside radius r, and the pressure of neutrinos, respectively. In equation H2-10() . the 
relativistic correction of magnetic pressure is ignored. Since the energy equation 
is formulated in the original code so as to obtain the internal energy density by 
subtracting the variations of kinetic energy and gravitational energy from the total 
work done by the matter pressure, the magnetic contribution to the kinetic energy 
must be also subtracted (magnetic pressure does no work on internal energy). Thus, 
the equation for energy conservation is replaced by 

dt ' dt\2 r 



dm r"^ dt\r 

P2 dt^ ' r ^ r 



-2T:e~^ — {P + Pu) - -rUq + -ATrrrF^ - tQ 



(2-13) 



where e is the specific internal energy, and Fi, and Q are the energy flux and the 
net heating rate of neutrinos, respectively. The general relativistic corrections from 
magnetic pressure are also neglected in the present study for simplicity. They are 
small enough in any case. 

Although the above modifications are sufficient in principle to calculate the 
evolutions of the magnetized wind, we have made other technical modifications. 
This is because of the numerical difficulties caused by adding the magnetic field 
terms to the original equations. When equations (12-10(1 and (|2-13() are employed to 
calculate the internal energy, numerical errors are accumulated, leading to unphysical 
solutions (such as decreasing entropy and negative internal energy). To avoid this 
problem, we have made two additional modifications. First, we solve the entropy 
equation with the heating and cooling processes of neutrinos taken into account, 

e-'^T^ = -muTQ. (2-14) 
at 

Here is the atomic mass unit. Second, instead of employing the energy density 
derived from equation (|2-13() , we use the entropy given by equation (|2-14() to obtain 
the energy density from the equation of state (EOS). It is true that equation (|2-14l) 
is not valid when a shock wave appears, but this is not a concern in this paper, since 
we investigate only the subsonic winds. In fact, the shock does not appear in this 
study. 

Apart from the modifications listed above, the code is identical with the one used 
by Ref.*^ The only exception is the boundary condition for the matter pressure, 
which are shown below for detail. We use the relativistic EOS table which was 
developed by Shen et al.^'*^ As the wind expands, the density becomes lower than 
the minimum value of the original EOS table. We have extended the table for lower 
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densities {p < 10^ g/cm^). Our treatment of neutrinos is simple. Following Ref.,*^ 
we assume that the neutrino distribution is constant in time and expressed as 

fuAr) = fu,AEu,-{Eu,)), (2-15) 

where fi,.^ is given by (the neutrino luminosity L^^. , the radius of the neutrino sphere 
R^^, and the average neutrino energy {Ey^)) as 

27r(l - x)Ly^ 

J^ir = /F \sp2 ' (2-16) 




(2-17) 

I - I 

The average energies are taken as follows, 

{EyJ = WMeV, (2-18) 

{EpJ = 20MeV, (2-19) 

{EyJ = 30MeV. (2-20) 

In the present study, we do not distinguish ^, r neutrinos and their anti-neutrinos 
with each other, so the identical average energy is adopted for them. The luminosity 
is assumed to be the same for all flavors and is taken as 

= lO^^ergs/s. (2-21) 

We assume that the radius of neutrino spheres R^. are all the same with that of 
PNS, R. 

In the present paper, we set the radius R = 10km and gravitational mass M = 
IAMq, which are typical values for observed NS. A very thin surface layer (initially 
~ 400 m) of PNS which contains a baryon mass of ~ 1O~^M0 is constructed on 
the numerical grid. Within this layer, 100 zones are uniformly spaced in the baryon 
mass coordinates (~ 10~^Mq per zone). 

Since the wind is expected to be blown from a nearly hydrostatic PNS surface, 
the initial density profile of matter is determined by solving the Oppenheimer-Volkoff 
equation in this surface layer by setting the density at the innermost mesh as 2 x 
lO^'^g/cm^ and assuming a constant value for the temperature T and the electron 
fraction in the whole layer. In the present study, we take T = 3MeV and Yg. = 0.25. 
By setting the matter pressure at the outer boundary (see equation ()2-25|) ) by hand, 
there is a little deviation from a perfectly static configuration. We run the code for 
a short time without heating and cooling processes of neutrinos until a complete 
hydrostatic state is obtained for the given outer boundary condition. We adopt 
the resulting configuration as the initial condition. The hydrostatic distributions 
obtained by this procedure are shown in figure ^ The electron fraction remains 
constant (=0.25) because no neutrino reactions are included. 
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We have assumed a homogenous configuration for the initial magnetic fields. 
Thus, even after adding them to the matter distribution obtained above, the initial 
condition remains hydrostatic. As mentioned previously, the direction of the mag- 
netic field is assumed to be perpendicular to the radial direction (Bg{m) = const, 
Br^Bfj, = 0). We adopt a magnetic field strength ranging from ~ 10^^ G (typical 
value for an ordinary NS) to ~ 5 x 10^^ G (suitable for a magnetar). Together 
with the non-magnetic model (model N), the models in this study are listed in Ta- 
ble ni The names of these models are characterized by their initial plasma beta f3 
= P/{B'^ /Stt) at the outer boundary. 



Model 


Bg{G) 




Pout 




Pb(dyn/cm'^) 


N 











1.0x10^^ 


lO^'M 


1.0x10^^ 


1.6x10® 


1.0x10^ 




1.0x10^^ 


lO^M 


1.0x10^^ 


1.6x10* 


1.0x10^ 




1.0x10^^ 


lOM 


1.0x10" 


1.6x10^ 


1.0x10 




1.Qx1Q22 


IM 


3.2x10" 


1.6x10 


1.0 




9.9x1Q22 




1. 0x10^^5 


1.6 


1.0x10" 


1 


9.0x10^^ 




1.4x10^'^ 


7.9x10"^ 


5.0x10" 


2 


8.0x10^^ 




2.2x10^^5 


3.2x10"^ 


2.0x10" 


2 


6.3x10^^ 




3.2x10^^^ 


1.6x10"^ 


1.0x10" 


2 


4.7x10^^ 


2DUM 


4.5x10^'^ 


7.9x10^2 


5.0x10" 


3 


2.6x10^^ 



Table I. Initial models. Column 1 gives the names of the models. Column 2 presents the strength 
of the initial magnetic field. Columns 3 and 4 display, respectively, the initial plasma betas in 
the innermost and outermost mesh points. In column 5, matter pressure Pb which appears in 
equation l|2-25|l (the outer boundary condition) is shown. Pb is set so as to give the same total 
boundary pressure (P + S^/Svr). 



We have employed free inner and outer boundary conditions except for radius 
and velocity. At the inner boundary, the velocity is set to be zero, and radius is fixed 
in time. 

Magnetic fields Bq are defined at the mesh interfaces in the numerical code. 
Their inner boundary condition is given as 



Bei = Bgi^i, 



(2-22) 



where the subscript I represents the innermost mesh interface. The outer boundary 
condition for magnetic fields is given as 



B, 



ON 



Bdex if Boex > and B^ex < Bqn-i 
BeN-i if Boex > BgN^i 
if Beex < 



(2-23) 



where the subscript N stands for the outermost interface, and B^x is the linear 
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extrapolation of magnetic field given as 

BoN-2irN-i - tn) + BoN-iirN - rN-2) „ ,s 

rN-1 - rN-2 

The third condition in equation H2-23p is imposed to avoid the negative magnetic field, 
although this situation has not occurred in our calculations. The outer boundary 
condition for the matter pressure is expressed in a similar way as 



Pn 



Pex if Pex > Pb and Pex < Pn~i 

Pn-1 if Pex > Pb and Pex > Pn-1 , (2-25) 

Pb if Pex < Pb 



where Pex is the linear extrapolation of matter pressure multiplied by a reduction 
factor 0.94, and is given as 

Pex = 0.94 X Pn-Arn-r^)+PnirN-rn-i) _ (2-26) 

'^n ^n— 1 

The subscripts n and N denote, respectively, the outermost mesh center and inter- 
face. Pf, is constant which varies among the models (see Table QJ). These values are 
chosen so as to give an asymptotic temperature suitable for r-process (T ~ O.lMeV) 
in the final phase of the wind evolution. In fact, we have set Pf, to give the identical 
total pressure {P + B'^ /8tt = lO^^dyn/cm^) at t = 50s. In order to suppress oscilla- 
tions at the boundary, we have limited the boundary pressures to give only negative 
gradients. The extrapolation has been employed instead unlike the previous study 
of Ref.^) so as to avoid an unrealistic accelerations at the outermost mesh point 
by the strong magnetic field [B > 10^^), which may produce a shock wave. The 
factor 0.94 which appears in equation (|2-26() was employed only just it minimized 
the deviation from the hydrostatic state when making initial conditions, and there 
are no physical backgrounds. It is noted that the prescription has little infiuence on 
the overall dynamics. 



§3. Numerical results 



We have computed all the models in Tabled by the same procedure given in §2. 
The only difference in each model is the magnetic field strength and the value of Ph, 
the boundary value of matter pressure. First, we show the numerical result of the 
non-magnetic wind case which was already studied by Ref.^^ (with a little change in 
the boundary conditions and numerical treatments of the basic equations). We have 
done this to check for the validity of our numerical method and to use the result for 
comparison with the magnetic wind solutions. 

Figures show the evolution of the non- magnetic model (model N) . The r- 
process in the neutrino-driven wind proceeds as follows. At the surface of the PNS, 
matter is initially composed of relativistic particles, free nucleons, and photons. 
The nascent PNS contracts and releases ~ lO^^ergs of gravitational energy mostly 
through neutrino production and emission over its Kelvin-Helmholtz cooling time 
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('^ 10 — 20seconds), and a small fraction of the energy is deposited at the surface 
layer via following reactions. 

fe + n^e~+p (3-1) 

+ P e"*" + n (3'2) 

Ui + e~ ^ Ui + e~ (3-3) 

Ui + e+ + e+ (3-4) 

i/i + e~ ^ Ui + e~ (3-5) 

z^i + e+ ;7j + e"^ (3-6) 

Ui + Ui ^ e~ + e"*" (3-7) 

Here the subscript i denotes the flavor of neutrino {i = e,^,T). Most efficient re- 
actions to deposit energy are the neutrino absorption on free nucleons f (|3-l|) and 
H3-2|) ). These reactions also determine the evolution of the electron fraction. When 
sufficient energy is deposited to overcome the gravitational attraction, an outflow 
(neutrino-driven wind) is generated. The wind is heated until it reaches ~ 50km. 
As the wind expands, the temperature decreases and when it drops below ~ IMeV, 
the free nucleons start to combine with each other and form a-particles. At this 
stage, neutrino reactions have ceased, and the subsequent evolution becomes adia- 
batic. When the temperature drops below ~ 0.5MeV, a-process occurs (e.g. triple 
a reactions), and seed nuclei of the r-process are formed until the temperature drops 
to ~ 0.2MeV. 

There are three key quantities in the wind scenario which determine the neutron- 
to-seed ratio. Those are the electron fraction (Ye) and entropy per baryon (s), and the 
dynamical time scale {T^iyn). The dynamical timescale is defined to be the duration 
of the a-process. For a robust r-process, lower values are favored for Ye and r^y^ 
because it implies that neutrons are rich and the duration of the production of 
seed nuclei is short. Higher entropy is favored for the following reason. When the 
system is still in high temperature regime (above ~ 0.5MeV), and nuclear statistical 
equilibrium (NSE) is maintained, the high entropy shifts the equilibrium to the state 
in which matter is composed mainly of a— particles and free neutrons with only a 
small fraction of nuclei. However, when the temperature becomes < O.SMeV, the 
system falls out of NSE. The high entropy implies a low density. Therefore, the 
a-process proceeds slowly, and most of the matter remains as a-particles. 

In model N, the wind expands rapidly ~0.1 s after the beginning of neutrino 
heating processes. As the temperature decreases, the a-processes occur, which corre- 
sponds to the location from '^50km to ~150km. The state is reached to the asymp- 
totic one with a constant temperature ~0.1MeV at ~0.2s. Rapid increases of the 
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electron fraction and entropy are seen when the wind is expanding from the surface 
of the PNS to ~50km. Note that a shght increase of entropy in the late phase of the 
wind evolution is caused by a breakdown of approximation for the heating by elec- 
tron and positron scatterings. This occurs because the equation we have used for the 
heating due to electron and positron scatterings is valid only when the temperature 
is high and, therefore, overestimates the increase of the entropy in the low temper- 
ature regime (see Ref.^-* for detail). This takes place in all the models commonly. 
However, this affects little the dynamics of the wind in this phase, and does not 
make a significant change in the dynamical time scale, electron fraction and entropy 
before the a-process stage. Thus, we can safely discuss the possibility of r-process 
by ignoring this false increase in entropy. Figure |31 shows the baryon mass density 
p, matter pressure P, and entropy per baryon s as a function of baryon mass. The 
density and pressure show significant decrease when the wind is expanding rapidly. 
Then the asymptotic state appears where p ~ lO^g/cm^ and P r-^ lO^^dyn/cm^ in 
the whole region. These features are identical with the previous studies done by 
Ref.8) 

The trajectories of models lO^M, ^^M, and are shown in figure ID As the 

initial magnetic field becomes greater, more rapid expansion can be seen. However, 
the later evolution shows only little deviation from model N (see figure [TJ, and the 
overall evolution for all models is nearly identical with that for the non-magnetic 
model. Figures [S] and IHl present, respectively, the evolutions of temperature and 
electron fraction for models lO^M, jqM, and 250^' ™ model N, the efficient 
decrease in temperature and increase in electron fraction take place when the matter 
is rapidly expanding. The values of temperatures and electron fractions at the same 
radius almost coincide among all the models including model N. This implies that the 
models with stronger magnetic fields show an earlier drop in temperature and rise 
in the electron fraction, and the later behavior becomes identical. The asymptotic 
values of Yg and T do not vary so much among all models (Ye ~0.45, T ~ O.lMeV). 

The distributions of the baryon mass density p, entropy per baryon s, matter 
pressure P as well as magnetic pressure B'^/Stt for these three models are given 
(as a function of the baryon mass) in figures ISl llUI The evolutions of baryon mass 
density and pressure are essentially same in all models (including model N). The 
density and pressure decrease as the wind expands and reach an asymptotic value, 
which is nearly identical among all models {p ~ lO^g/cm^, p ~ lO^^dyn/cm^). The 
evolutions of the magnetic pressure and entropy are also similar among the models 
except for the asymptotic value. As the wind expands, the initially homogeneous 
magnetic pressure decreases and reaches an asymptotic value. Unlike the density and 
pressure, the asymptotic value varies largely in the region. The final distribution of 
magnetic pressure has a positive gradient, and the outer region is ~2 orders of 
magnitude larger than the inner region. Since the neutrino heating processes occur 
up to r -^SOkm, the entropy increases until the wind reaches that radius. The models 
with initially strong magnetic fields tend to have lower asymptotic values compared 
to the models with weaker fields. This tendency can be observed particularly in 
the outer region of the wind. As mentioned above, the more rapid evolutions of 
these values {p, P, s, and B'^/Stt) seen in the models with stronger magnetic fields 
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are results of the more rapid expansion. These behaviors wih be discussed in more 
detail in the next section. 

The three key quantities (s, Tdyn-, and Yf>) for the r-process are summarized for all 
the models in Table HII Since the outer and inner regions are affected by the boundary 
conditions, the average values shown in the table are obtained by discarding the 
inner and outer 5 mesh points. The entropy in the table is the average value at 
the beginning of the a-process (T~0.5MeV). The false increase of entropy caused 
by the breakdown of the approximation in the late phase is ignored. The dynamical 
time scale Tdyn is calculated as the duration, in which the temperature decreases 
from O.SMeV to 0.2MeV, for each trajectory and averaged. As noted above, we can 
not see any significant difference in the trajectories of the magnetic wind cases and 
non-magnetic wind case. Even in the models with magnetar-scale field strengths, 
only a slight change comes out (the trajectories of model 200^ together with model 
N are displayed in figure [71). Therefore, it is no surprise that the evolution of the 
variables such as temperature, electron fraction, density, pressure and entropy are 
also nearly identical (see figures 01151 and ISllUf) . As a result, even strong magnetic 
fields give remarkably little change in the three key parameters (see Table 111}) . We 
will discuss the reason for this in the next section in greater detail. 



Model 


s(fef,) 








N 


137 


4.63x10" 


-2 


0.453 


lO^'M 


137 


4.57x10" 


-2 


0.453 


10^ M 


136 


4.57x10" 


-2 


0.453 


lOM 


136 


4.55x10" 


-2 


0.453 


IM 


137 


4.55x10" 


-2 


0.453 




136 


4.41x10" 


-2 


0.453 




135 


5.27x10" 


■ 2 


0.453 




133 


5.35x10" 


-2 


0.453 




129 


4.01x10" 


-2 


0.453 


2ro^ 


128 


5.45x10" 


-2 


0.453 



Table 11. Three key quantities for each model. These values are the average over the whole wind 
after discarding the inner and outer 5 mesh points. 

§4. Discussions 

Our results show that the homogenous magnetic field has only a little effect on 
the wind dynamics. In particular, the three key quantities for the r-process do not 
vary much from those for the non-magnetic wind. This is true even for the very 
strong magnetic fields. In this section, we will analyze this behavior in more detail. 

The reason why the dynamics is little affected by magnetic fields is understood 
from the evolutions of magnetic fields. Figures Iimi2l show the evolutions of the 
plasma beta, which is defined to be the ratio of matter pressure to magnetic pressure. 
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for models lO^M and 2So^" these models, the plasma beta increases rapidly 
from the initial values and the magnetic pressure becomes negligible for the wind 
dynamics. This tendency can be explained analytically as follows. 

From equation (|2-8|) . the variation rate of the magnetic pressure is given as 

1 dB"^ V dX , 

2- -2-17, (4-1) 



5| dt r dt' 

where d/dt denotes the Lagrange derivative, and w is a radial fluid velocity which is 
defined as 

-I- 

Ignoring the neutrino heating, which is important only at the earliest stage, we can 
evaluate the variation rate of matter pressure approximately as 

IdP I dp 

P dt p dt 

where 7 is the adiabatic index. Ignoring further general relativistic corrections, e'^ 
can be approximated as 

^ ( dm \ ~^ 1 



^ , , dX 
From the above equation, — is written as 

dt 

dX I dp V 

37 ~ — -J^-^-- (4-6) 
dt p dt r 

With the use of the continuity equation and equation i)4-6|l , equations (|4-1|1 and H4-3|) 
become, respectively 

(4-7) 



1 dBl 








= 2^-2( 




Bl dt 


r \ 




1 dP 


V 


/ dv 


pIF ~ 


-27 7 

r 


V dr 



(4-8) 

If we locally approximate the velocity profile by a power-law, v (x r", the variation 
rates are evaluated as 

J^f|l,_p + 2„K-S (4.9) 

1 dP 

p— cx-(2 + a)7r"-i. (4-10) 

Since the adiabatic index 7 is approximately 4/3, we finally obtain the criterion that 
if a < 1, the plasma beta is decreased. 
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Figure El shows the radial profile of velocity for model N when the wind has 
settled to a quasi-steady state. Though this is a snapshot for this model, all the 
models have a similar configuration. If we apply the above criterion in this figure, we 
find that the region in which the plasma beta should decrease is the outer region (r > 
60km). It is noted that the assumptions we have made to derive the above condition, 
that is, the adiabaticity and nonrelativisity of the flows, are clearly satisfied in the 
region. Thus, after the wind expands beyond ~60km, the plasma beta will decrease. 
Since we are considering the region exterior to the PNS surface, the Newtonian 
gravity is not so bad an approximation in all region, and equation ()4-9() is valid even 
in the inner region (r < 60km). On the other hand, equation H4-10() is not appropriate 
in this inner region, since the heating cannot be ignored and the adiabatic condition 
is not satisfied. However, it is obvious that the heating processes tend to increase the 
matter pressure and, as a result, increase the plasma beta more than in the adiabatic 
flow. Thus, although the obtained criterion is not valid, the plasma beta is expected 
to increase in the inner region where a > 1. 

From the above discussions, the wind is expected to expand initially with in- 
creasing plasma beta up to r ~60km, and thereafter with decreasing plasma beta. 
Our numerical results actually agree with the expectation. This can be confirmed 
by seeing figures [TT1I12I Since the increase of the plasma beta occurs more rapidly 
than the following decrease, the wind tends to be matter-pressure-dominant in the 
entire evolution. Thus, even models with very strong magnetic fields give only a lit- 
tle effect on the wind evolution. In the following, we will further analyze, although 
small, effects of magnetic fields on the dynamics, focusing on model 200^' 

At the beginning of expansion, the magnetic pressure is dominant, as can be 
seen from figures ISl llOl and its negative gradient pushes the matter outward. This 
causes an early rapid expansion seen in figure [Tj Although the magnetic pressure 
accelerates matter for a while, the effects of the magnetic pressure soon becomes 
significantly small (since the plasma beta decreases rapidly). Then, as they approach 
the asymptotic state, the gradient of the magnetic field changes signs and starts 
to decelerate matter (see figures 111112(1 . Note that although magnetic pressure is 
small at this stage, it still has some leverage on the dynamics to cancel the initial 
accelerations. This causes the later catch-up of non-magnetic wind model N seen in 
figure [71 

The subsequent distribution of the magnetic pressure is a direct consequence of 
the initial assumption on the magnetic field. The initial homogeneity of magnetic 
fields implies that the matter near the surface has a greater magnetic fiux per unit 
mass compared with the matter deeper inside, since the density is higher there. 
As the wind evolves to an asymptotic state with a constant density in radius, the 
frozen-in magnetic field has a larger magnetic fiux densities in the outer region. In 
another word, the final magnetic field configuration has an inverted profile of the 
initial matter density distribution. 

This can be understood more precisely from equation ()2-8() and figure ^ By 
using the same approximation that we have made previously, we can approximate 
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equation (|2-8|) as 



Be{ti,m) 



r{ti,m) p{ti,m) 
r{to,m) p{to,m) 



Be{to,m), 



(4-11) 



for arbitrary to and ti. Taking to and ti, respectively, as the initial time and the 
time at which the asymptotic state is reached, we find that r(to,?7i), r(ti,m), and 
p{ti,m) depend on m very weakly. Bg{tQ,m) is also identical on all mesh points. 
Thus, the ratio of magnetic pressures between two mesh points in the asymptotic 
state can be approximated as 



From the initial density distribution shown in figure ^ and the evolutions of density 
and magnetic pressure shown in figures ISllUl this relation is verified. 

As a result of the initial expansion and the later deceleration of the flow men- 
tioned above, the key quantities for the r-process except for the electron fraction show 
only a little change among the models. On the other hand, this change depends on 
the initial location of mass element in the atmosphere, which can be roughly divided 
into (i) the outer region, (ii) the central region and (iii) the inner region as follows: 

(i) The outer region (1.2 x IQ-'^Mq < m < 1.6 x IQ-^Mq) 

Since the outer region has the lowest plasma beta, the early rapid expansion occurs 
most efficiently. Consequently, as matter passes through the heating region rapidly, 
the resulting entropy s becomes lower compared with the non-magnetic wind. This 
effect can be seen by comparing figureOlwith figures ISl llUI Moreover, the deceleration 
takes place earlier on in the outer regions because the density gradient is initially 
larger than the inner regions. This occurs during the a-process stage, and makes the 
dynamical time scale Tdyn longer. 

(ii) The central region (4.0 x W'^Mq < m < 1.2 x W^Mq) 

Though their plasma beta is not so low as in the outer region, the initial rapid ex- 
pansion takes place also in this region and the entropy s becomes lower than for 
the non-magnetic wind. Unlike the outer region, however, the deceleration occurs 
much later after the wind has already passed through the a-process region (~50km 
to ~150km). Thus, the dynamical time scale becomes shorter. The three key quan- 
tities in the central region are displayed for all the models in Table ITTll in which we 
have taken average by excluding the outer and inner 25 mesh points. The reduction 
of the dynamical time scale in the models with strong magnetic fields can be ob- 
served in the table. However, the reduction is very small and, moreover, the entropy 
decreases at the same time. Thus, even if we focus only on this region, a successful 
r-process will not take place. 

(iii) The inner region (m < 4.0 x IO^^Mq) 

The inner region also experiences the acceleration and deceleration stages. However, 
because of their large initial plasma beta, their effects are too small and the key 
quantities for the r-process are essentially unaffected. 

As mentioned previously in §2, our calculations approximate the flow near the 
equator with initial magnetic fields perpendicular to the radial direction. Our results 



Bg{ti,m2y p{to,miy 
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Model 


s{h) 








N 


130 


3.37x10" 


'2 


0.453 




130 


3.37x10" 


-2 


0.453 


10^ M 


130 


3.37x10" 


-2 


0.453 


iUM 


1 on 
ioO 


■i.ol X iU 


-2 


U.4o3 


IM 


130 


3.35x10" 


-2 


0.453 




129 


3.36x10" 


■2 


0.453 




129 


3.33x10" 


■2 


0.453 




133 


3.08x10" 


■2 


0.453 




124 


2.83x10" 


'2 


0.453 




121 


2.78x10" 


-2 


0.453 



Table III. The three key quantities for each model. The quoted values are the average over the 
region excluding the inner and outer 25 mesh points. 



can be applied to the places in the wind where the deviation from the radial motion 
is small, that is, regions where the radial force dominates over the transverse force. 
The angular range that satisfies this condition can be estimated as follows. For an 
arbitrary initial field-configuration, if the transverse velocity is small compared with 
the radial velocity, the magnetic fields are approximately determined by equations 
H2-7|) - ()2-9() . In the nonrelativistic limit, these equations become 

Br{ti,m) = !^ij^S,(to,m), (4-13) 



Bg{ti,m)= f /' : Beito,m), (4-14) 

r{to,m)p{to,m) 

r{ti,m)p{ti,m) 

B^{ti,m) = — — -B^{to,m). (4-15) 

r{to,m)p{to,m) 

The evolution of the magnetic fields is determined by the initial magnetic fields, 
radius, and density, as well as by the present radius and density. The typical values 
for the initial radius and density, and the final radius and density in our models are 
10km, lO^'^gcm"^, 500km, and 2xl0'^gcm~^, respectively. If the initial magnetic field 
is uniform and perpendicular to the the equatorial plane, the radial force dominates 
over the transverse force roughly in the angular range between 65° ^ < 115°. This 
occupies a significant fraction of the total solid angle, and our results are applicable 
to a rather large region below and above the equatorial plane. 



§5. Summary 



In this paper we have investigated the effects of magnetic fields on the neutrino- 
driven wind and explored the possibility as an r-process site by performing numerical 
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simulations. The homogenous distribution is assumed for the initial magnetic fields, 
and ten models with a wide range of field-strength have been computed. 

The results do not give the dynamical behavior expected by Ref.^^^ This may 
be the consequence of the assumed wind configuration, which has a tendency to 
increase the plasma beta rapidly and the matter-pressure is dominated in most of 
the wind evolutions. Hence, the magnetic pressure has little chance to influence the 
dynamics. 

The evolutions of weak magnetic field models are identical with non-magnetic 
model. The models with very strong magnetic fields {B > lO^^G) show some, though 
not large, variation. For about 0.5 x 10"^ Mq of wind, the reduction of the dynamical 
time scale has been found. However, this change is very small and accompanied by 
the decreases of entropy. Thus, unfortunately, the situation favored for a successful 
r-process is not realized in our models, even with magnetar-scale field-strengths. 

Since our study has been limited to initially homogeneous magnetic fields, other 
distributions of magnetic field should be investigated. Note, however, that the initial 
surface structure will be changed substantially then, since the magnetic pressure 
will determine the hydrostatic force balance. Multi-dimensional MHD simulations 
of magnetar formation may shed a light on this issue. It is also self-evident that 
more realistic and multidimensional investigations will be needed further, since the 
field-tension that is not considered in this paper may play an important role. Even 
in that case, very strong magnetic fields > lO^^G will be needed to influence the 
wind dynamics. Last but not least, magnetic rcconnections might be interesting, 
since the dissipation of magnetic fields will heat the matter, accelerating the wind. 
However, this is much beyond the scope of this paper. 
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Fig. 1. The initial configuration of the surface layer. Radius r, baryon mass density p, matter 
pressure p, temperature T, electron fraction Ve, and entropy per baryon s are given as a function 
of baryon mass coordinate. 
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Fig. 3. The baryon mass density p, pressure P, and entropy s in model N as a function of baryon 
mass. Each line shows the distribution at t=2.53 x lO'^s, 1.00 x 10"^s, 1.22 x 10"^s, 1.58 x 10"^s, 
1.80 X lO^^s, 2.12 X 10~^s, 2.48 x lO^^s, 3.26 x 10"^s, and 5.00 x 10s. In the innermost region 
in the above figures, p and P decrease and s increases monotonically. The distribution of s at 
t=5.00 X 10s is not displayed, since it shows a false increase (see the text for detail). 
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Fig. 6. Evolution of the electron fraction for model lO^M (top left panel), (top right panel), 

and (bottom panel). Every 5 mesh points are plotted. 
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Fig. 7. Trajectories for model 25olV[(thick solid lines) and trajectories for model N(thin solid lines). 
Every 20 mesh points are plotted. 
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Fig. 8. The baryon mass density p, pressure P, entropy s, and magnetic pressure B'^/Stv for model 
lO^M as a function of baryon mass. Each line shows the distribution at t=2.53 x 10~*s, 1.00 x 
IQ-^s, 1.22 X lO'^s, 1.58 x 10~^s, 1.80 x 10"^s, 2.12 x 10~^s, 2.48 x 10"^s, 3.26 x 10"^s, and 
5.00 x 10s. In the innermost region in the above figures, p, P, and B^/Stt decrease and s increases 
monotonically. The distribution of s at t=5.00 x 10s is not displayed, since it shows a later false 
increase (see the text for detail). 
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Fig. 9. The baryon mass density p, pressure P, entropy s, and magnetic pressure B'^/Stv for model 
as a function of baryon mass. Each line shows the distribution at t=4.26 x 10~^s, 1.02 x 
10"^s, 1.22 X 10"^s, 1.56 X 10"^s, 1.85 x 10"^s, 2.16 x 10~^s, 2.53 x 10"^s, 3.40 x 10~^s, 
and 5.00 x 10s. In the innermost region in the above figures, p, P, and B^/Stt decrease and s 
increases monotonically. The distribution of s at t=5.00 x 10s is not displayed, since it shows a 
false increase (see the text for detail). 
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Fig. 10. The baryon mass density p, pressure P, pressure s, and magnetic pressure B'^/Stt for model 
2^M as a function of baryon mass mesh. Each line shows the distribution at t=7.24 x 10~^s, 
1.00 X 10~^s, 1.22 X 10~^s, 1.56 x 10"^s, 1.83 x 10"^s, 2.13 x 10"^s, 2.34 x 10"^s, 2.84 x 10~^s, 
and 5.00 x 10s. In the innermost region in the above figures, p, P, and /8tv decrease and s 
increases monotonically. The distribution at t=5.00 x 10s is not displayed, since it shows a false 
increase (see the text for detail). 
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Fig. 11. Evolution of the plasma beta for model lO^M. Each line shows the distribution at t=2.53 x 
10-'*s, 1.49 X IQ-^s, 4.12 x lO'^s, 1.00 x lO'^s, 1.22 x lO'^s, 1.58 x lO'^s, 2.12 x lO'^s, 
2.33xl0~^s, and 5.00x10s. In the innermost region in the above figure, /3 increases monotonically 
in time. 
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Fig. 12. Evolution of the plasma beta for model ^qM.- Each line shows the distribution at t=7.24 x 
lO'^'s, 1.41 x 10"^s, 4.10 X 10"^s, 1.00 x 10"^s, 1.22 x 10"^s, 1.56 x lO^^s, 1.83 x 10"^s, 
2.13x 10~^s, and 5.00x 10s. In the innermost region in the above figure, P increases monotonically 
in time. 
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Fig. 13. Velocity profile for model N at t = 2.33 x 10"^s. 



